Thermoelectric Response Near the Density Driven Mott Transition 
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We investigate the thermoelectric response of correlated electron systems near the density driven 
Mott transition using the dynamical mean field theory. 



Thermoelectric effects which are responsible for the di- 
rect conversion of heat energy into electrical energy and 
vice versa have recently received renewed attention [Q. 
So far optimal thermoelectric response has been achieved 
in traditional doped semiconductors but recent improve- 
ments in both the synthesis of correlated electron systems 
and the theoretical methods for treating them have gen- 
erated new interest in the thermoelectric properties of 
correlated materials. 

The thermoelectric performance of a material, reflects 
its ability to convert applied voltages into tempera- 
ture gradients while minimizing the irreversible effects 
of Joule heating and thermal conduction. This can be 
quantified by the so called dimensionless figure of merit 
denoted here by ZtT: 



Field Theory (DMFT) USEffl, which has successfully 
described many aspects of the physics of three dimen- 
sional transition metal oxides. The goal of our study is 
to understand qualitatively the thermoelectric response 
of a doped Mott insulator and to derive explicit formulas 
for the transport coefficients which are valid at very low 
and high temperatures. Using these insights and numer- 
ical results we discuss the figure of merit of this partic- 
ular kind of materials. An early numerical calculation 
of the Seebeck coefficient in the large d Hubbard model 
appeared in ref. . 

We consider perfect periodic solids described by the 
TV-fold degenerate Hubbard model: 
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Here a is the electrical conductivity, S the thermopower 
(or Seebeck coefficient) and T the temperature, k is the 
electronic contribution to the thermal conductivity and 
Kr, is the lattice thermal conductivity. In weakly corre- 
lated systems the Seebeck coefficient can be interpreted 
as the logarithmic derivative of the conductivity with re- 
spect to the Fermi energy. In strongly correlated electron 
systems S is a more difficult quantity to interpret, and is 
this the subject of this letter. 

We consider a system of fcrmions doped away from a 
Mott insulating state, where the magnetic correlations 
are weak so that the magnetism is not the driving force 
behind the metal to insulator transition. This situa- 
tion is realized experimentally in titanate and vanadate 
perovskite compounds [Q. For example, Lai^Y^TiOa 
is ferromagnetic for x near 1 and anti-ferromagnetic at 
small values of x but in all these compounds the Neel and 
Curie temperatures are quite small, of the order of 130 
K and less || so we focus on the paramagnetic phase of 
the doped Mott insulator. 

The experimental motivation for our study was mainly 
provided by the work of Y. Tokura's group which has 
demonstrated that in the class of ternary perovskites, one 
can control to a large extent various parameters such as 
orbital degeneracy, bandwidth, and carrier concentration 
to provide an experimental realization of the filling driven 
Mott transition 

To approach this problem we use the Dynamical Mean 



and ignore electron-phonon interactions. The index a 
can be thought of as a spin or an orbital index, and will 
run from 1 to N. The single band case corresponds to 
N = 2. We now summarize the relevant aspects of the 
DMFT. The transport properties are obtained from a 
Green's function with a frequency dependent but mo- 
mentum independent self energy: 



G(k,u) =G(e k ,u) = 
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where e& is the dispersion relation. The self-energy, is 
computed by solving an Anderson impurity model in a 
bath described by a hybridization function A(iu). Re- 
garded as a functional of the hybridization function it 
obeys the self consistency condition |]|| ■ 



dsG(s,iuj)D(e) (4) 
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where D(e) is the bare density of states, D(e) = J^k 3{ £ ~ 
€k). The transport coefficients that govern the electrical 
and thermal responses of the model are given in terms of 
current-current correlation functions. Within the DMFT 
they reduce to averages over the spectral density p(e, lo) 
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In this expression the relevant information about the bare 
band structure is contained in the spectral density, p, and 
the transport function <&, 



(7) 



The DMFT expressions bear superficial similarity with 
recent work of Mahan and Sofo |||, if we identify their 
transport function which we denote S M s(£) as: 
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(8) 



Ref. M stressed the relevance of correlated materials 
to thermo-electricity and suggested that the optimal re- 
sponse is obtained when E MS takes the form of a delta 
function like peak at the Fermi level. The DMFT allows 
us to derive explicit expressions for E MS starting from 
microscopic Hamiltonians. Near the Mott transition, the 
low energy spectral function does assume a delta like form 
at low temperatures. We find however, that the optimum 
of the figure of merit is achieved at high temperatures, 
when the quasi-dclta-function-likc resonance at the Fermi 
level is absent. 

In electronic systems near a Mott transition there are 
two widely separated energy scales: the bare bandwidth 
D, which sets the scale for incoherent excitations and 
the coherence temperature T co k where Fermi-liquid-like 
properties begin to be observed. T co -^ vanishes as we 
approach half filling. As a result when we decrease the 
temperature, starting from very high temperatures we 
expect two crossovers to take place, one when fc^T « D 
and the other when T w ^oh- We assume in this analysis 
that the interaction energy U is much larger than the bare 
bandwidth, i.e. we are close to the the U = oo limit. 

At very low temperature the Fermi liquid picture is 
applicable and momentum space offers a natural descrip- 
tion of the transport processes. The only states that 
contribute come from a narrow region around the Fermi 
surface, which in the local mean-field picture shows up 
as a narrow Kondo like resonance near the Fermi level. 
At very high temperatures the entire Brillouin zone is 
equally important and a real-space picture of the trans- 
port is more appropriate. The transport is then entirely 
due to incoherent motion of charge carriers. We discuss 
these two asymptotic regimes below. 

Low temperature regime: In this regime we use 
Fermi liquid ideas to parameterize the transport coeffi- 
cients in terms of a few parameters, which contain all the 
effects of the interactions: 7 = Jm£(0) is the scattering 



rate at the Fermi level, Z = (1- 



-■^f^L^or 1 is the quasi- 
particle residue at the Fermi-level and p — p — ReT,(0) 
is the effective chemical potential. Close to the Mott 
transition there is only one effective scale controlling the 
low-tempcraturc physics and therefore the scattering rate 
behaves as A(k B T) 2 / Z 2 D where A is a dimensionless 
constant |L0| ]. Performing a low-temperature expansion 
of equation (ra) and substituting into (m, we obtain: 
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We find that the thermopower and the figure of merit do 
not depend on 7 and are given by 
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Here the numbers E n are given by, 
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Numerical calculation gives Eq = 0.82 and E% = 1.75. 
Note that the Wiedeman-Franz law holds here as it does 
for transport dominated by impurity scattering. How- 
ever, since the scattering rate is energy dependent around 
the Fermi-surface the ratio of aT and k is not equal to the 
classical Lorenz number. The physical content of equa- 
tion (jll]) is transparent. Correlations enhance the figure 
of merit relative to that of a non-interacting system with 
the same density of states by a factor of i which in 
this context can be thought of as the square of the mass 
enhancement. This factor is expected to be large and 
in fact diverges as we approach the density driven Mott 
transition in Lai-^Sr^TiOa 0. Note however that ex- 
pression (|llj) is only valid in the low-temperature regime 
because it is restricted to (3 ZD 3> 1. Thus the figure of 
merit will be very low in the low-temperature regime un- 
less the logarithmic derivative of the transport functions 
becomes appreciable. Unlike the density of states how- 
ever, the transport function does not have any Van Hove 
singularities and the only singular points in the logarith- 
mic derivative are at the band edges where the transport 
function vanishes. 

High temperature regime: To describe this region 
we observe that the spectral function shifted by the value 
of the chemical potential 
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converges to a well-defined, shape centered around w = 
as the temperature tends to infinity. This agrees with a 
rigid-band interpretation of the lower Hubbard band, ex- 
cept that the carriers in this band are completely incoher- 
ent near the band edge. The high-temperature behavior 
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of the chemical potential can also be found analytically 



and is given by: (3/i = ln(- 



We can now obtain 



^ iV(l-n) ' 

the leading high-temperature behavior of the transport 
coefficients by inserting the scaling form ([l3|) with p tem- 
perature independent and expanding to lowest order in 
(3 the equations for the coefficients A n . The results are 
parameterized in terms of the moments 
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which we evaluate numerically. To leading order in [3 we 
find: 
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Near the Mott transition the density and degeneracy de- 
pendence of the moments 70 and 72 is given by: 7 m = 
7 m (l — n + jf) 2 where the 7 m 's are constants, 70 ~ 0.05 
and 72 « 0.01. The factor in parenthesis is simply the 
total integrated spectral weight of the Green's function. 
The high-temperature equation that we get for the ther- 
mopower corresponds to the well known Heikes formula 
jll| generalized for degeneracy N. Comparison with nu- 
merical solutions of the dynamical mean field equations, 
reveals that for the resistivity, the high temperature ex- 
pansion formula is quite accurate over a very wide tem- 
perature range and breaks down only at temperatures 
of the order of This is surprising, since the high 

temperature expansion is a priori only valid for (3D <C 1. 
This result gives some insight into the origin of the lin- 
ear resistivity which was observed in early studies of the 
Hubbard model in infinite dimensions p2| . Furthermore 
if one keeps the next term in the high-temperature expan- 
sion of the thermopower, one can fit the numerical data 
over a comparable region. Notice that the thermopower 
close to the Mott transition is hole- like which agrees with 
the picture of holes in a paramagnetic spin background 

& 

The high-temperature expression for the figure of merit 
is to lowest order in (3 given by: 
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where nn — ^P-. Here the lattice contribution to the 

^ an 

thermal conductivity has been included since the elec- 
tronic contribution tends to zero in the high tempera- 
ture region. As temperature increases the figure of merit 
increases monotonically to a constant value that goes lin- 
early with the bandwidth, D, of the system. At any finite 



temperature however the equation above gives a figure of 
merit that increases monotonically with D to a maxi- 
mum which is at a bandwidth larger than ksT and thus 
outside the region of validity of our formula. Thus we 
conclude that the optimum figure of merit is obtained 
when the bandwidth is of the order of the temperature. 

We now turn to quantitative calculations of the ther- 
mopower in the intermediate temperature range which is 
the range most pertinent to experiments and possible ap- 
plications. The best characterized system, in the class of 
materials that we seek to describe is the Lai-^Sr^TiOa 
system with x small. We model this system with a 
Hubbard model on a three dimensional hypercubic lat- 
tice with half bandwidth D — 0.5 eV and interaction 
strength U = 2.0 eV in the Mott insulating end of the 
series (x = 0). We take into account the x dependence of 
the bandwidth by using the fact that the Ti-O-Ti bond- 
angle, 9, changes with doping. The bandwidth then de- 
pends on 9 through D{6) = D(180°) cos 2 9. To select 
the bond-angles corresponding to a given doping we use 
data for ab-plane bond-angles from G3| . The results from 
the calculations of S, using IPT, with this choice of pa- 
rameters is displayed in Fig. [l]. We notice here that in 
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FIG. 1. Thermopower vs. temperature for several values 
of doping in the Lai-^Sr^TiOa system. 

the high temperature end of the data the values of the 
thermopower are actually close to what Heikes formula 
predicts thus indicating that the high temperature ther- 
mopower can be used as a rough estimate of the carrier 
density. At lower temperature we see a marked decrease 
in the magnitude of the thermopower which heralds the 
transition into the Fermi-liquid regime where the ther- 
mopower is small and electron-like. We note however 
that for doping values lower than 0.1 the thermopower 
is already hole-like above 300K which indicates that the 
quasiparticle peak has almost disappeared and that the 
carriers are completely incoherent. This regime does not 
appear for smaller values of the interaction U when the 
system is below the Mott transition. The thermopower 
near the Mott transition therefore seems to be a sensitive 
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indicator of the character of the carriers (i.e. localized vs. 
itinerant). 

Since bandwidth and carrier concentration, can be con- 
trolled to a large degree in this class of materials we use 
our results to assess the prospects for a large figure of 
merit in this class of systems. We are not aware of any 
measurements of the thermal conductivity in the titanate 
compounds so for the purpose of obtaining an order of 
magnitude estimate of the figure of merit we use a value 
of 2.0W/mK for kl based on measurements carried out 
recently in the manganite oxides Jl5[|, 

The results from the numerical calculation of the fig- 
ure of merit in the intermediate temperature region are 
displayed in Fig. |. The data for T = 0.05 and T = 0.10 
were obtained with iterated perturbation theory (IPT) 
modified for finite doping jl6). For the T = 0.20 data we 
use the infinite U non crossing approximation (NCA). 
This approach has also been shown to give results in 
good agreement with exact methods at high tempera- 
tures |0,Q8|. The behavior of the figure of merit can 
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FIG. 2. Thermopower vs. temperature for several values 
of doping. Calculated with IPT using U = 4.0. 

be understood qualitatively as follows. It vanishes at a 
density no(T) where the Seebeck coefficient vanishes. At 
very high-temperature this occurs at rig = ^ N accord- 




ing to equation (|16|). As temperature is lowered no(T) 
moves towards half filling, and is given by the condition 
T con (7i) « T . This is due to the formation of the coher- 
ent quasi-particles which give a negative contribution to 
the thermopower. For densities lower than tiq(T) the sys- 
tem is Fermi-liquid-like and we find a low figure of merit, 
closer to the Mott transition at temperature higher than 
T con (n) the figure of merit is larger. 

The lowest temperature results show little density de- 
pendence below no- This is because both the logarithmic 
derivative of the transport function and the quasiparti- 
cle residue depend essentially linearly on the doping and 
therefore the figure of merit varies little with density. The 
T = 0.1 data seems to be displaying a similar trend as 



the T = 0.05 data but the system is essentially out of the 
Fermi-liquid regime at that temperature. 

To conclude, we have investigated the thermoelectric 
coefficient and the thermoelectric figure of merit near the 
density driven Mott transition. We provided simple ex- 
pressions for the transport coefficients of this model, ob- 
taining a qualitative understanding of the thermoelectric 
coefficient. In the light of our results, a large figure of 
merit in this class of systems seems rather unlikely. It 
would require very small doping, and a very small value 
of the bare half-bandwidth parameter D, to be able to 
access the high temperature regime T > D where the 
figure of merit is of order unity. We have argued how- 
ever, that the thermopower is a sensitive probe of the 
degree of itineracy of the carriers, and an experimental 
investigation of the thermoelectric response in the region 
where large mass enhancement has been observed [|| is 
highly desirable. 
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